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Abstract 

We find a large family of solutions to the Dirac equation on a manifold of Gi holonomy 
asymptotic to a cone over S* 3 x S 3 , including all radial solutions. The behaviour of these 
solutions is studied as the manifold developes a conical singularity. None of the solutions 
found are both localised and square integrable at the origin. This result is consistent with 
the absence of chiral fermions in M-theory on the conifold over S 3 x S 3 . The approach here 
is complementary to previous analyses using dualities and anomaly cancellation. 
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1 Introduction 



Eleven dimensional supergravity cannot give a chiral four dimensional theory upon compacti- 
fication on a smooth manifold One way to resolve this discrepancy with experiment is by 
compactifying on a manifold with boundary (2|, ||. Recently, work has centred on obtaining chiral 
fermions by the alternative approach of compactifying on singular manifolds. As supergravity is 
no longer valid near the singular points, the arguments make essential use of M-theory. Given 
that M-theory remains unformulated, this means they work via dualities with string theories. 
Approaches include duality with intersecting D6-branes and orientifolds of type IIA string theory 
0, ^, P and fibrewise duality with Heterotic string theory compactified on certain Calabi-Yau 
manifolds Q . The duality results are corroborated by Witten's analysis of anomaly cancellation 
at conical singularities ||. 

A recurrent theme in these arguments is compactification on manifolds of G2 holonomy (|| 
for mathematical background, [|n]] for a pedagogical introduction). This is unsurprising as su- 
persymmetry plays an important role in duality arguments and G2 holonomy is the condition on 
the internal seven dimensional manifold for compactification from eleven to four dimensions to 
preserve N = 1 supersymmetry. The preservation of N = 1 supersymmetry is phenomenolog- 
ically appealing and G2 manifolds were considered in early work on Kaluza-Klein supergravity 
1 11, [13]. The possibility of further obtaining chiral theories with non-abelian gauge symmetry 
via M-theory is even more phenomenologically interesting. 

Parallel to the renewed interest in the physics of G2 holonomy compactifications, there has 
been a lot of recent work on constructing and classifying G2 metrics. All the examples known 
are cohomogeneity-one metrics on noncompact spaces that are asymptotic to generalised cones 
|l|, U, |l|, Ul, 0, H], EI EH H H- The metrics constructed are generalisations of the metics 
found a decade ago by p3|, |24]]. They are asymptotic to cones over SU(3)/U(1) x U(l), CP 3 
and S 13 x S 3 . No explicit metrics are known on compact G2 manifolds because these cannot 
have continuous isometries except for trivial S 1 factors. The appearance of chiral fermions is 
a phenomenon localised near the conical singularity and so considering noncompact manifolds 
should not be problematic. 

The original argument against obtaining chiral theories from eleven dimensions 0] is cir- 
cumvented by the fact that the non-abelian gauge groups under which the chiral fermions are 
charged do not come from isometries of the internal manifold, but from massless branes wrap- 
ping collapsed homology cycles in the internal manifold |2f|, |26|| . It would be nice to see these 
new massless fermionic degrees of freedom emerge as homology cycles collapse directly within an 
eleven dimensional framework, without using dualities. Further, for the G2 manifold asymptotic 
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to a cone over S 3 x S 3 , call it X, the situation regarding chiral fermions is not yet understood, 
although other aspects of M-theory on this space have been studied in some detail (e.g. |27|). 
There are various reasons for this: (a) The anomaly analysis of ||] is not applicable as for this 
space H 2 (X;U(l)) = 0. (b) The G 2 manifolds considered in 0] were M 3 bundles over a self- 
dual four- fold whilst X is an R 4 bundle over S 3 . (c) Whilst the cones over SU(3)/U(1) x U{\) 
and CP 3 were interpreted in Q as the M-theory lift of intersecting D6-branes, the analogous 
interpretation for the S 3 x S 3 case is not clear. 

The line of attack in this work will be to look for solutions to the massless Dirac equation on 
manifolds on G2 holonomy and see how the behaviour of the solutions changes as the manifold 
developes a conical singularity. We use the Dirac equation and not the Rarita-Schwinger spin-3/2 
equation because the degrees of freedom we are interested in do not originate from supergravity 
but from wrapped branes. It turns out that the G2 manifold on which the Dirac equation is most 
tractable, allowing us to find all radial solutions and several non-radial solutions, is precisely the 
one asymptotic to a cone over S 3 x S 3 . Thus the information we gain through this approach is 
nicely complementary to the information already obtained through dualities. Solving the Dirac 
equation on other G2 manifolds could also be interesting, but is not done here. 

In §2 we write down the massless Dirac equation on the G2 manifold and find all radial 
solutions by directly solving the equation. In §3 we find some non-radial solutions by relating 
solutions of the massless Dirac equation to solutions of the massless Klein-Gordon equation via 
the covariantly constant spinor on the G2 manifold. The qualitative behaviour of the solutions 
does not change in the singular limit and none of the solutions have the property of being both 
localised at the origin and square integrable at the origin, properties which we should expect 
for degrees of freedom corresponding to the chiral fermions. In §4 we discuss these results in 
the context of chiral fermions appearing in singular M-theory compactifications. We suggest 
that our results are consistent with the absence of chiral fermions arising in the singular limit of 
M-theory compactified on a G2 holonomy manifold asymptotic to a cone over S 3 x S 3 . 

2 Dirac equation and all radial solutions 

We will work with the simplest G2 holonomy metric asymptotic to a cone over S 3 x S 3 which 



has metric 




ds 2 = a(r) 2 dr 2 + /3(r) V - -£*) 2 + 7 (r) 2 S i S i , (1) 
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where the a % and the are left invariant one- forms on the two copies of 5 3 , that is da 1 = 
— ie ljfe cr-? A a k and dTi 1 = —\e % i k YJ> A £ fe . The radial functions are 

*(rf = (1 - jf)- 1 ; /?(r) 2 = ±r 2 (l - jf); 7 (r) 2 = lr 2 . (2) 

The range of r is ro < r < oo. It is useful to consider the vielbeins 

e° = adr- i = (5{a l - -E*); & = 7 E\ (3) 

We will use tangent space indices i,j,k = 1,2,3 and i,j,k = 1,2,3 = 4, 5,6 and also a,b,c... = 
0...6. Spacetime indices will be [i,v.... 

The massless Dirac equation corresponding to the metric (Q) can be read off from the calcu- 
lation of the connection one-form in |24|] , 

T a D a ^ = T a (e a - - A J C a Y b }j V> 



+^ fer M^ " + * = o, (4) 



•8 7 *j 167 2 « ' 87" 

where as usual T ao = | [r a , 1^] and e a are vector fields dual to the vielbeins @, < e b , e a >= 5 b a . 
We use a convention for the gamma matrics such that the Clifford algebra is {r a ,rfe} = —25 ao - 
To solve this equation we need to choose a representation for the Euclidean Clifford algebra 
in seven dimensions. For this algebra there is a representation in which all the matrics are real. 



It is given by [ESL 29] 



(r a )st = C as t + 8 as dt7 - bat^sT, (5) 

where the spinor indices s,t = 0...7 and c as t are zero if s = 7 or t = 7 and otherwise are the 
totally antisymmetric octonion structure constants, with 

C016 = 0352 = 0)43 = Cl42 = C135 = C236 = C456 =1. (6) 

In this representation, the Majorana or reality condition is particularly simple 

r = v>- (7) 

The Dirac operator @) takes Majorana (real) spinors to Majorana spinors. Therefore by ex- 
pressing an arbitrary spinor in terms of two Majorana spinors 

t/j = -ifjW +iV {2) , (8) 

we may consider the Dirac equation for ip^ and ip^ separately. 
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For the remainder of this section we take tp to be Majorana. We will thus obtain half of the 
solutions. However, the other half will be of exactly the same form but with a pure imaginary 
rather than real coefficient. Writing out ifi in components 



■0 — (hi hi hi hi hi hi hi h 



(9) 



where fi-.-h are rea l functions. The Dirac equation thus becomes eight coupled equations for 
eight real functions. These equations are only tractable if one makes the radial ansatz 



eiif) = e-ip = 0, 

for i = 1, 2, 3. That is, we restrict the spinor to depend on r only. 
Substituting (||) into the Dirac equation (||) one obtains the ei, 



(10) 
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(11) 



Note that these equations form three groups; for {/ 2 , h}, for {/ 4 , f 5 } and for {fi, / 3 , / 6 , /§}. 
To discuss the solutions, we need to distinguish the cases r 7^ and r = 0. In each case we 
obtain eight linearly independent solutions which, when combined with the other eight solutions 
of the same form but pure imaginary coefficient, are all the radial solutions to the Dirac equation 
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2.1 Solutions with r ^ 

Define f = r — tq. We are interested in whether the solutions are square integrable, i.e. L 2 , near 
the origin. For this we need the determinant of the metric (|l]) which near the origin is J~g ~ f 
as f — > 0. The behaviour of the eight solutions is summarised in the following table. 



Independent solns. 


As r — ► 


As r — > oo 


L 2 at origin? 


Localised? 


{/2,,/V}; {/4,/s}; fz, h, fs} 




~ f 


Yes 


No 


{/2,jV}; {.A, is}; {h, f-i, h, fs} 


~ f- 3 / 2 




No 


Yes 


{/lj /3j /6, /8} 






Yes 


No 


{/l, /3> /6, /8} 


~ f- 3 / 2 




No 


Yes 



Where the notation in the table means that, for instance, there is a solution in which f± 
and /s go as a constant near the origin and as r at infinity and all other radial functions are 
zero. The second and third columns show only the asymptotic power-law behaviour, not the 
numerical coefficients, which in general are not the same for the nonzero radial functions of 
each solution. The general radial solution is a linear superposition of these eight solutions with 
complex coeffients. None of the solutions are both 1? at the origin and localised at the origin. We 
are not interested in the behaviour at infinity, although it is easy to see that none of the solutions 
are L 2 everywhere, in agreement with Lichnerowicz's theorem for noncompact manifolds. 

2.2 Solutions with r = 

The form of the determinant at the origin is now ^fg ~ r 6 as r — > 0. The following table 
summarises the behaviour of the eight solutions. 



Independent solns. 


As r — > 


As r — > oo 


L 2 at origin? 


Localised? 


{/2, /r}; {/4,/5}; {/l>/3,/6,/8} 


r 


~ r 


Yes 


No 


{h, h}; {/4,/s}; h, fa, fs} 


-7 




No 


Yes 


{fl, /3) hi /8} 


~ r° 




Yes 


No 


{/l> /3; h, h} 


~ r -6 


~ r~ 6 


No 


Yes 



The power-law asymptotics are in fact exact solutions in this case. We see that the asymptotic 
behaviour at infinity did not depend on ro, which is obvious from (Q) anyway. Whilst the 
behaviour at the origin does change, it is qualitatively the same and still none of the solutions 
are both L 2 at the origin and localised at the origin. The possibility of the Lichnerowicz theorem 
not holding existed in this case because the manifold is singular. 
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3 Non-radial solutions 



Some non-radial solutions may be obtained by using relationships between eigenfuntions of dif- 
ferent differential operators. This is possible because of the existence of covariantly constant 
spinors on special holonomy manifolds |J. The idea is developed in |3(]] for four dimensional 
hyperKahler metrics and in p3] for massive eigenfunctions on G2 and Spin (7) metrics. The mass- 
less case we need is slighty more subtle. We will find solutions for the massles Dirac equation 
from solutions to the massless scalar equation. Suppose that 4> satisfies 

A0 = -V V> = 0, (12) 

and that 77 is the covariantly constant spinor, Dr/ = 0. Then 

if (V a (/>)r a ?7 = then ip = 4>t] solves r a V a ^ = 0, 

if (V ^>)r ry / then rp = (V a c/>)r a 7/ solves Y a V a ^ = 0. (13) 
One can also obtain solutions of the Dirac equation from solutions of the massless vector equation 



1 24], but this is no easier to solve than the Dirac equation itself and so is not useful to us here. 
The first thing we need to calculate is the covariantly constant spinor. This can be done 



|24, 31] by imposing the projection conditions 

(r 04 - r 23 ) r? = (r 05 - r 3 i) v = (r 06 - r 12 ) v = o, (u) 

and then requiring 77 to have constant norm rfr/ = 1 and to be Major ana as in (0). This first 
condition is consistent because D(rfr]) = and the reality condition is necessary because the 
supersymmetry generator must be real for degrees of freedom to match up in the supergravity 
theory. All this determines 77 up to a sign, which can be chosen. One obtains 

7^ = ^(1,0,-1,0,0,1,0,1). (15) 

So the covariantly constant spinor has constant components. This will be generically true for 
cohomogeneity one metrics. The immediate consequence is that all r dependence of ip will come 
from the scalar solution 4>. By substitution it is easy to check that Dr\ = as required (cf. the 
constant solutions of §2.1 and §2.2 above). 

The massless Klein-Gordon equation on the manifold can be written 

= - (eZety/gdvt) = 0, (16) 

where ea are the inverse vielbeins. This expression for the Hodge-de Rham operator is well suited 
for exploiting the symmetries of cohomogeneity-one metrics ([!]). It has been used to study the 



6 



Schrodinger equation on the Atiya-Hitchen metric in |}^] and on the Eguchi-Hanson metric in 
1 33]. For the metric ([[]), the equation becomes 



' -a r (r 6 (l - r o 3 /r 3 ) 2 r 0) + ^J 4 J,</> + 4 ! + ^M<f> = 0, (17) 



r 6 (l — rjj/r 3 ) r 2 4 * r r 2 1 — r|]/r 3 r 2 ' 

where Jj are the vector fields dual to X* and the J[ are dual to a 1 . It follows that these vectors 
satisfy the algebra [Jj, Jj] = e^^Jk and [JJ, Jj] = HjkJ'k- This means that as differential operators, 
JiJi, J'iJ'i and Ji J[ are simultaneously diagonalisable. Further, their eigenvalues are well known 
from elementary angular momentum theory. Acting on the simultaneous eigenfunctions <3?™ n we 
have 

J i J i *T = -rn(m + l)*% n ; A^T = ~<n + 1)*™ ; J^T = st^T , (18) 

where m, n G t^Z + U {0} and {s, t} & {m, n} + Z with — m < s < m and — n < t < n. Note that 
the half integers are allowed because the group in question is SU(2) and not SO^). We are not 
interested in the precise functional form of the angular eigenfunctions. Details and references 
can be found in [^]. One now decomposes the function into modes of the form <f> = $(r)$^ n . 
The equation (TH]) becomes 



d 2 ^ / ^ i 6r 2 d ^ in ( m(m + 1) + st i i iN 4r 3 — r$ 



This equation has four regular singular points, at r 3 = rjj and at infinity. It does in fact have 
a solution in terms of hypergeometric functions, but these are only valid for r < rg whilst the 
range of r is r > vq. 

This reduction to a single ordinary differential equation was only possible because of the 
relatively large symmetry of the metric (|]). More general G2 metrics with less symmetry would 
have resulted in a Klein-Gordon equation in which the differential operators were not all simul- 
taneously diagonalisable. 

The behaviour at infinity is independent of rg and we have that as r — > 00 

$( r ) ^ r ^(5±v / 25+48Q) ) (2Q) 

where Q = m(m + 1) + n(n + 1) + st > 0. In the radial case, Q = 0, the two possible behaviours 
are just the | ip |~ r° and | ip |~ r -6 solutions which we found in the previous section. For the 
second solution we use | if? |~ V r <j) ~ d r (j) as discussed above. The case of Q / will be discussed 
below. The behaviour at the origin depends on whether ro is nonzero. 
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3.1 Solutions with r ^ 

The behaviour as r — > for the non radial modes depends on whether the integer n is nonzero. 
If n / then as f — » 

^.^K^)). (2i) 

The power of r is nonzero and therefore | ij) |~ r l l 2 df4> ~ f 2 ( 2=t v /l+4 ' 1 ("'+ 1 )) _ ^11 Q f these 
solutions either are not localised or are not square integrable at the origin. Recall that the 
metric goes as yjg ~ f. The effect of having angular dependence on these solutions is to increase 
the degree of divergence at the origin of the localised solutions and increase the degree to which 
the non-localised solutions increase away from the origin. 
And if n = then as f — > 

$(f) ~ k 1 f° + k 2 f- 1 . (22) 

These correspond to spinors behaving like | ip |~ f° and | ip |~ f~ 3 / 2 . The first does not decay 
and the second is not I? at the origin. This set of solutions does not change its radial behaviour 
as the angular dependence changes. The behaviour found is again consistent with the results of 
the previous section. 

3.2 Solutions with r = 

The solution as r — > is the same as in the r — > oo limit of (|2CJ|). We now consider the effect of 
angular dependence. All of the solutions are nonconstant in r and so 

| V |~ 9 r ~ r ^-(7± v / 25+48Q) i (23 ) 

Again, all of these solutions either are not localised or are not square integrable at the origin, 
because Q > 3/4. The metric now goes as ^fg ~ r 6 . Again, the effect of including angular 
dependence is to make the solutions divergent at the origin more divergent and the solutions 
that grow away from the origin grow faster. 

4 Discussion 

We have found solutions that are constant, solutions decaying away from the origin and solutions 
growing away from the origin. None of these satisified the two properties of being both localised 
and square integrable at the origin. The first of these is expected for degrees of freedom emerging 
as a homology cycle collapses and the second is necessary for the solution to make sense quantum 
mechanically. Allowing angular depedence worsens the situation in all the solutions we found. 
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The results admit three possible interpretations. Firstly, that the degrees of freedom corre- 
sponding to chiral fermions localised near the singularity correspond to nonradial solutions of 
the Dirac equation that we have not found. Secondly, that solving the Dirac equation is not the 
way to exhibit these degrees of freedom in any case, or that a modified Dirac equation including 
gauge connections is required. Thirdly, that the G2 manifold asymptotic to a cone over S 3 x S* 3 
does not give chiral fermions in the singular limit, unlike the cones over SU (3)/C7 (1) x U (1) and 
CP 3 . 

The third of these possibilities is the most interesting. It links nicely with the results from 
dualities and anomaly cancellation discussed in the introduction. It also fits with the picture 
of extra massless modes coming from collapsed membranes and the fact that the G2 metric 
asymptotic to a cone over S 3 x S 3 has zero second homology group. Solving the Dirac equation 
on other G2 metrics that are known to give chiral fermions appears to be more challenging, but 
could lend further support to this possibility. Another interesting question is to consider the 
effect of including gauge field backgrounds in the Dirac equation. 



Many thanks to Gary Gibbons for useful suggestions. Thanks also to Ruben Portugues, 
James Sparks and Tibra Ali for discussions. The author is funded by the Sims scholarship. 
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